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Abstract 

The Synchronized Molecular-Dynamics simulation which was recently pro¬ 
posed by authors [Phys. Rev. X 4, 041011 (2014)] is applied to the analysis 
of polymer lubrication between parallel plates. The rheological properties, 
conformational change of polymer chains, and temperature rise due to the 
viscous heating are investigated with changing the values of thermal conduc¬ 
tivity of the polymeric liquid. It is found that at a small applied shear stress 
on the plate, the temperature of polymeric liquid only slightly increases in 
inverse proportion to the thermal conductivity and the apparent viscosity of 
polymeric liquid is not much affected by changing the thermal conductivity. 
However, at a large shear stress, the transitional behaviors of the polymeric 
liquid occur due to the interplay of the shear deformation and viscous heat¬ 
ing by changing the thermal conductivity. This transition is characterized by 
the Nahme-Griffith number Na which is defined as the ratio of the viscous 
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heating to the thermal conduction at a characteristic temperature. When the 
Nahme-Griffith number exceeds the unity, the temperature of polymeric liq¬ 
uid increases rapidly and the apparent viscosity also exponentially decreases 
as the thermal conductivity decreases. The conformation of polymer chains 
is stretched and aligned by the shear flow for Na < 1, but the coherent 
structure becomes disturbed by the thermal motion of molecules for Na > 1. 
Keywords: Multiscale modeling, Polymeric liquid, Lubrication, Rheology 


1. Introduction 

Molecular dynamics (MD) simulations are used to investigate the materi¬ 
als properties and internal structures of the complex liquids and soft matters. 
In MD simulations, the dynamics of all the model molecules which compose of 
the materials are calculated. The material properties and internal structures 
of the materials are analyzed from the results of the MD simulations. The 
MD simulation is applicable for a wide variety of complex liquids and soft 
matters consist of the complicated molecules such as the polymer molecules 
and the biological molecules. However, due to the enormous computational 
time to calculate the dynamics of all the molecules, the MD simulation is not 
yet applicable to problems which concern mesoscopic and macroscopic scales 
motions far beyond the molecular size, where the dynamics of molecules are 
strongly coupled to the macroscopic transports of the complex liquids. At 
microscopic scales, the macroscopic transports are usually ignored and the 
MD simulation is performed under a certain ideal environment. To develop 
a simulation technology to analyze the molecular dynamics coupled with the 
macroscopic transports for the complex liquids is a significant challenge from 
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both a scientific and engineering point of view. Multiscale modeling is a 
promising candidate to address this challenge. 

The multiscale simulation for the flow behaviors of complex fluids was 
first advanced in the CONNFFESSIT approach for polymeric liquids by Laso 
and Ottinger 4, [5|, 6], where the local stress in the fluid solver is calculated 
using the microscopic simulation instead of using any constitutive relations. 
The GENERIC approach is also presented for the non-isothermal polymeric 
flows by making important corrections and clarifications to the CONNFFES¬ 


SIT scheme. [7] The strategy exploited in the CONNFFESSIT approach is 
also introduced into the heterogeneous multiscale modeling (HMM), which 
was proposed by E and Enquist,^] as a general methodology for the ef¬ 
ficient numerical computation of problems with multiscale characteristics. 
HMM has been applied to various problems, such as the simple polymeric 
flow 9], coarsening of a bina ry p olymer blend 10] and the channel flow of a 
simple Lennard-Jones liquid ll[. The equation-free multiscale computation 
proposed by Kevrekids et al. is also based on a similar idea and has been 


applied to various problems. 
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3 , [13] De et al. proposed the scale-bridging 
method, which can correctly reproduce the memory effect of a polymeric 
liquid, and demonstrated the non-linear viscoelastic behavior of a polymeric 
liquid in slab and cylindrical geometries. 14j, [15] The multiscale simulation 
for polymeric flows with the advection of memory in two and three dimen¬ 
sions was developed by Murashima and Taniguchi.flG. 17, 181 We have also 
developed a multiscale simulation coupling of MD and transport equations. 
The multiscale method was first developed for simple fluids 
quently extended to polymeric liquids with the memory effect 
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Recently, we proposed the synchronized molecular dynamics (SMD) method 
in which the multiscale method was extended to treat the coupled heat and 


momentum transfer of complex liquids. 


24] 


In the previous study, the SMD method was applied to the polymer lu¬ 
brication generating the viscous heating between parallel plates. The rheo¬ 
logical properties, conformation of polymer chains, and temperature rise due 
to the viscous heating were investigated for various applied shear stresses on 
the plate. We found that an interesting transitional behavior of the confor¬ 
mation of polymer chains occurs with a rapid temperature rise for a large 
applied shear stress due to the coupling of the shear deformation and heat 
generation under the shear flow. 

In this study, we investigate the effects of changing the thermal con¬ 
ductivity of the polymeric liquid on the behaviors of rheological properties, 
conformation of polymer chains, and heat generation in the polymer lubrica¬ 
tion between parallel plates. In the following, we first describe the problem 
considered in the present paper. The simulation method is briefly explained 
after the presentation of the problem. The SMD simulation of polymer lu¬ 
brication is performed, and the results are discussed; these results are mainly 
the rheological properties and the coupling of the conformations of the poly¬ 
mer chains with heat generation under shear flows. Finally, a short summary 
is given. 


2. Problem and Method 

We consider a polymeric liquid contained in a gap of width H between 
parallel plates with a constant temperature T 0 (see Fig. QJa)). The polymeric 
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A constant shear stress is applied to the upper plate. 



MD cells are asigned to each small fluid element to calculate the 
local quantities. 
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(a) Geometry of problem 


(b) Time evolution 


Figure 1: The schematic of the geometry of the problem (a) and the time evolution (b). 
This figure is reprinted from Ref. n 


liquid is composed of short Kremer-Grest chains 25] of ten beads, in which 


all of the bead particles interact with a truncated Lennard-Jones potential 
defined by 

. 6 " 


Ulj (r) = 


4e 

0, 


(d 12 -(d‘ 


+ e i 


(r < 2 1 / 6 ct), 
(r > 2 1 / 6 ct), 


( 1 ) 


and consecutive beads on each chain are connected by an anharmonic spring 
potential, 

U F (r) = ~k c F% In 

with /c c =30e/cr 2 and i? 0 =1.5cr. The polymeric liquid is in a quiescent state 
with a uniform density po and a uniform temperature To before a time t — 0. 
Hereafter, the y -axis is perpendicular to the parallel plates, and the bound¬ 
aries between the upper and lower plates and the polymeric liquid are located 
at y = H and 0, respectively. The upper plate starts to move in the x direc¬ 
tion with a constant shear stress c 0 at a time f=0, while the lower plate is 
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at rest. 

The macroscopic behavior of the polymeric liquid is described by the 
following transport equations: 


dv x 


de 

po dt 


d(Jxy 

d 2 T 

a, ' f< ~ 


(3a) 

(3b) 


where v a is the velocity, a a p is the stress tensor, e is the internal energy 
per unit mass, and 7 is the shear rate, i.e., 7 = dv x /dy. Hereafter, the 
subscripts a, /3, and 7 represent the index in Cartesian coordinates, i.e., 
{a,/ 3 , 7 }E{x,y,z}. Here, we assume that the macroscopic quantities are uni¬ 
form in the x and z directions, d/dx=d/dz= 0 , and the density of the poly¬ 
meric liquid is constant. Fourier’s law for a heat flux with a constant and 
uniform thermal conductivity A is also considered in Eq. (13b|h Note that the 
therrna 
general 


conductivity of polymeric liquids is anisotropic under shear flows in 


26l |27l ]28l 129], and some experimental studies have reported that the 
linear stress-thermal relation between the stress tensor and thermal conduc¬ 
tivity holds. 0 , 0 , 0 , Q However, in the present study, we only consider 
the isotropic thermal conductivity as the first step because the effect of shear 
thinning of the viscosity is thought to be more crucial to viscous heating un¬ 
der strong shear flows than that of the anisotropy of the thermal conductivity. 
We also assume that the velocity and temperature of the polymeric liquid 
are the same as those of the plates at the boundaries, i.e., the non-slip and 
non-temperature-jump boundary conditions. 

The effect of viscous heating is estimated using the ratio of the first and 
second terms in Eq. (l3b|) to be cr 0 T H 2 /AAT 0 . Here, T is the gross shear 
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rate of the system, which is defined by the ratio of the velocity of the upper 
plate v w to the width of the gap H,T — v w /H , and AT 0 is a characteristic 
temperature rise for the polymeric liquid. In the present problem, we consider 
a characteristic temperature necessary to substantially change the viscosity of 
the polymeric liquid, i.e., AT 0 = \rj 0 /(dr} 0 /dT 0 )\, where r/ 0 is the characteristic 
viscosity of the polymeric liquid at a temperature of Tq. Thus, the Nahme- 
Grifhth number Na, defined as 


Na = 


a 0 TH 2 


(4) 


A|cllog(77o)/c)To|- 1 

represents the effect of viscous heating on the changes in the rheological 


properties. 


34j Usually, in lubrication systems and in high-speed process¬ 


ing operations with polymeric liquids, the Nahme-Griffith number is not 
negligibly small because of the large viscosity and the small thermal conduc¬ 
tivity of the polymeric liquid. 34] For example, when a lubrication oil in a 
gap with a width of 10 /jrn is subjected to shear deformation with a strain 
rate of 1 x 10 4 s _1 , the Nahme-Griffith number is estimated to be Na > 0.1. 
Thus, the rheological properties of the lubricant in such micro devices must 
be significantly affected not only by the large velocity gradient but also by the 
temperature increase caused by local viscous heating. To predict the rheo¬ 
logical properties of the polymeric liquid in these systems, one must consider 
the temperature variation in Eq. (l3b|) coupled with Eq. (l3alh 

The viscous heating is one of the most fundamental and important prob¬ 
lems in the lubrication systems. However, it is difficult to analyze the behav¬ 
iors of polymer chains under the coupling of strong shear flows and viscous 
heating with using the conventional simulation technologies. As is seen in 
Eq. [3 the effect of viscous heating on the temperature and behaviors of 
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polymer chains is enhanced as the system size H is large. Thus, the thermal 
lubrication is one of the typical problems that the macroscopic transports 
and microscopic molecular dynamics are strongly interacted each other. We 
solve this problem by using the Synchronized Molecular-Dynamics (SMD) 


simulation which was proposed recently by authors 


24] 


In this paper, we briefly explain the SMD method. The molecular dynam¬ 
ics cells are assigned to small fluid elements with a interval Ax and calculate 
the local stresses, temperatures, and conformations of polymer chains in each 
MD cell by using the non-equilibrium molecular dynamics (NEMD) simula¬ 
tions according to the local shear rates 7 . (see Fig. [D(b)). The NEMD 
simulations are performed independently in each MD cell for a certain time 
duration At, but at each time interval At, the NEMD simulations are syn¬ 
chronized to satisfy the macroscopic heat- and momentum-transport equa¬ 
tions Eq. (J3]). Thus, the couplings of the shear flow, viscous heading, and 
conformational dynamics of polymer chains are autonomously calculated in 
each MD cell with satisfying the macroscopic heat- and momentum-transfers. 


3. Results 

In this study, we investigate the effect of changing the thermal conduc¬ 
tivity A of the polymeric liquid on the rheological property, conformation of 
polymer chains, and temperature rise. We carry out the SMD simulations 
with changing the value of the thermal conductivity A. (In the previous 
study 24j, a fixed value of the thermal conductivity is only considered while 
the applied shear stress to the upper plate <To is varied variously.) Here¬ 
after, we measure the length, time, temperature and density in units of a, 




r 0 = a Jma 2 /e, e/fcg, and m/a 3 , respectively. Here, ks is the Boltzmann 
constant, and m is the mass of the LJ particle. In the following simulations, 
the density of the polymeric liquid is fixed to be po = 1, and the temperature 
of the plates and the width of the gap between the plates are fixed to be 
T 0 = 0.2 and H = 2500, respectively. At this density p 0 and this temper¬ 


ature To, the conformation of the bead particles 
and results in complicated rheological properties. 


Decomes severely jammed 
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351] Two values of the 


applied shear stress on the upper plate a 0 are considered in this study, i.e., 
ao=0.01 and 0.05. The SMD simulations are performed with A =50, 100, 
150, 200, and 300 for a 0 =0.01 and A =50, 75, 100, 125, 150, 200, and 500 
for a o = 0.05, respectively. In the following, we present quantities averaged 
for a long period of time at the steady state, in which the shear stress is 
spatially uniform and the time derivative of the local internal energy, i.e., 
the left-hand side of Eq. (j3b|) . is negligible. 

Figure [5] shows the dependency of the apparent viscosity p, defined as p = 
<7o/T, on the thermal conductivity A and the spatial average of temperature 
T. For a small applied shear stress, i.e., <7o = 0.01, the apparent viscosity 
only slightly decreases as the thermal conductivity A decreases (the inverse 
of thermal conductivity A -1 increases). The average temperature T increases 
linearly as the function of the inverse of thermal conductivity A -1 , but the 
variation is very small; the temperature rise is at most 2.5% of the wall 
temperature in the figure. The facts that the variations of apparent viscosity 
p and average temperature T are small and the average temperature T is a 
linear function of the inverse of thermal conductivity A -1 can be explained by 
the energy transport equation, Eq. (l3bl) . for small Nahme-Grifhth numbers; 
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Figure 2: The apparent viscosity rj , defined as rj = cto/T, and the spatial average of 
the temperature T as functions of the inverse of thermal conductivity A^ 1 for the applied 
shear stresses (a) do = 0.01 and (b) cr 0 = 0.01, respectively. In (b), inset shows the 
spatial average of temperature as the function of A -1 for large thermal conductivities. The 
horizontal dotted lines on the left vertical axises in both figures indicate the viscosities for 
the infinitely large thermal conductivity, A^ 1 —> 0. In (b), the vertical arrows show the 
points at which the Nahme-Griffith number Na equals unity. 
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when the Nahme-Griffith number is small the viscosity of the polymeric liquid 
is not much affected due to the viscous heating, so that the local viscosities 
and shear rates are uniform between the parallel plates and the temperature 
rise is inversely proportional to the thermal conductivity. In Fig. [2]Ja) , the 
Nahme-Griffith number is at most 0.06. 

For a large applied shear stress, i.e., a o = 0.05, the behaviors of the 
apparent viscosity and average temperature are quite different from those 
for a small applied shear stress. For <7o = 0.05, the Nahme-Griffith number 
of the polymeric liquid can be larger than unity, Na > 1, at small thermal 
conductivities. In Fig. [2]^b), the points at which the Nahme-Griffith number 
equals unity are indicated by the vertical arrows and the Nahme-Griffith 
number is larger than unity on the right side of the vertical arrow. It is 
seen that the average temperature rapidly increases when the Nahme-Griffith 
number exceeds the unity, and the apparent viscosity also rapidly decreases 
as the function of the inverse of thermal conductivity. When the Nahme- 
Griffith number exceeds the unity, the viscosity substantially deceases due 
to the viscous heating, so that the shear velocity increases and the viscous 
heating is enhanced. This positive feedback causes the rapid changes of the 
apparent viscosity and average temperature for small thermal conductivities. 

Figure [3] shows the conformational changes in polymer chains as a function 
of the inverse of thermal conductivity. Here, the bond-orientation tensor Q a p 
is defined as, 


jy X/ ]\r 

1v p , • 1 v b 

^ chain 


N h -l 

E 


1 , Enin Enin 

3 =1 


(5) 


where N p is the number of polymer chains in each MD cell, iV b is the number 
of bead particles in a polymer chain, bj for 1 < j < Ay — 1 is the bond vector 
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Figure 3: The spatial average of the bond-orientation tensor of the polymer chains Q a p 
as a function of the inverse of thermal conductivity A -1 for the applied shear stresses (a) 
cto = 0.01 and (b) cto = 0.01, respectively. In (b), the vertical arrows show the points at 
which the Nahme-Griffith number Na equals unity. 

between consecutive beads in the same chain, and b min is the distance at which 
the sum + Up(r) has a minimum and is calculated to be 6 m ; n — 0.97. 

As the polymer chains are stretched and aligned by the shear flows driven 
by the applied shear stress to the upper plate, the difference of the xx- 
and yy-components of the bond orientation tensor, Q xx — Q yy , and the xy 
component of the bond orientation tensor Q xy becomes large. It is seen that, 
for a small applied shear stress (Fig. [3](a)), the conformations of polymer 
chains are slightly stretched and aligned by the shear flows. The variations 
of each component of the bond orientation to the thermal conductivity are 
very small; the relative variations are at most 0.3% for Q XX)Vy and 5% for 

Qxy 

For a large applied shear stress (Fig. [3](b)), the transitional behaviors 
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occurs around the point at which the Nahme-Griffith number equals unity. 
For large thermal conductivities (on the left side of downward arrow), the 
conformations of polymer chains are highly stretched and aligned by the 
shear flows, ffowever, for small thermal conductivities (on the right side 
of downward arrow), the alignment of the polymer chains is disturbed due 
to the thermal motion of molecules and the conformation of the polymer 
chains becomes more isotropic although the shear velocities are enhanced as 
the temperatures increase. This indicates that the transition between the 
regime dominated by the shear flow and dominated by the thermal motion 
of molecules occurs by changing the thermal conductivity. 


4. Summary 


We carry out the SMD simulations for the thermal lubrication of a model 
polymeric liquid composed of short chains between parallel plates. The SMD 
simulations for the same problem were also carried out in Ref. 24], where 
the behaviors of the polymeric liquid for various applied shear stresses are 
investigated. In the previous study we found that the transitional behav¬ 
iors occur in the dynamics of polymer chains by changing the applied shear 
stress; when the Nahme-Griffith number is small, Na < 1, the conforma¬ 
tion of polymer chains is stretched and aligned by the shear flow, but when 
the Nahme-Griffith number exceeds unity, Na > 1, the coherent structure 
becomes disturbed by the thermal motion of molecules. 

In the present study, the effects of changing the thermal conductivity of 
the polymeric liquid on the rheological properties, conformation of polymer 
chains, and temperature rise due to the viscous heating are investigated. It 
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is found that at a small applied shear stress on the plate, i.e., cr 0 = 0.01, 
the temperature of polymeric liquid only slightly increases in inverse pro¬ 
portion to the thermal conductivity and the apparent viscosity of polymeric 
liquid is not much affected by changing the thermal conductivity since the 
Nahme-Grifhth number is very small. However, at a large shear stress, i.e., 
<Jo = 0.05, the Nahme-Grifhth number of the polymeric liquid can be larger 
than unity, Na > 1, when the thermal conductivity is sufficiently small. 
When the Nahme-Grifhth number exceeds unity, the temperature of poly¬ 
meric liquid increases rapidly and the apparent viscosity exponentially de¬ 
creases as the thermal conductivity decreases. The conformation of polymer 
chains is stretched and aligned by the shear how for Na < 1, but the co¬ 
herent structure becomes disturbed by the thermal motion of molecules for 
Na > 1. Thus, in this study we demonstrate that the transitional behaviors 
which were found in the previous study also occur by changing the thermal 
conductivity under a fixed applied shear stress. 
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